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Abstract 

The production of a stochastic background of relic gravitational waves 
is well known in various works in the literature, where, using the so called 
adiabatically-amplified zero-point fluctuations process it has been shown 
how the standard inflationary scenario for the early universe can in prin- 
ciple provide a distinctive spectrum of relic gravitational waves. In this 
paper, it is shown that a weak modification of General Relativity produces 
a third massive polarization of gravitational waves and the primordial 
production of this polarization is analysed adapting the adiabatically- 
amplified zero-point fiuctuations process at this case. 

The presence of the mass could also have important applications in 
cosmology as the fact that gravitational waves can have mass could give 
a contribution to the dark matter of the Universe. 

At the end of the paper an upper bound for these relic gravitational 
waves, which arises from the WMAP constrains, is also released. 



1 Introduction 

Recently, the data analysis of interferometric gravitational waves (GWs) de- 
tectors has been started (for the current status of GWs interferometers see 
[BEllSllHISlllllTllH]) and the scientific community aims in a first direct detec- 
tion of GWs in next years. 

Detectors for GWs will be important for a better knowledge of the Universe 
and also to confirm or ruling out the physical consistency of General Relativity 
or of any other theory of gravitation [9l [lOl [TU [l2l [IBl [l^. This is because, in 
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the context of Extended Theories of Gravity, some differences between General 
Relativity and the others theories can be pointed out starting by the linearized 
theory of gravity [9l[l0l[l2l[Il]- In this picture, detectors for GWs are in principle 
sensitive also to a hypotetical scalar component of gravitational radiation, that 
appears in extended theories of gravity like scalar-tensor gravity and high order 
theories [I2l[lll[16l[l2l[18l[19l[20l[2ll|22], Brans-Dicke theory [22 and string 
theory [24] . 

A possible target of these experiments is the so called stochastic background 
of gravitational waves [25 t IMI^ [28 l [29 l [30] . 

The production of the primordial part of this stochastick background (relic 
GWs) is well known in the literature starting by the works of |251[26] and |271I28). 
that, using the so called adiabatically-amplified zero-point fluctuations process, 
have shown in two different ways how the standard inflationary scenario for the 
early universe can in principle provide a distinctive spectrum of relic gravita- 
tional waves, while in [291 ISD] the primordial production has been analyzed for 
the scalar component admitted from scalar-tensor gravity. In this paper, it is 
shown that a weak modiflcation of General Relativity generates a third massive 
polarization of gravitational waves and the primordial production of this polar- 
ization is analysed adapting the adiabatically-amplified zero-point fiuctuations 
process at this case. We have also to emphasize that, in a recent paper [40], such 
a process has been applied to the same theory which we are going to discuss 
in the present work. But, in [40] a different point of view has been considered. 
In that case, using a conform analysis, the authors discussed such a process in 
respect to the two standard polarizations which arises from standard General 
Relativity. In the present paper the analysis is focused to the third massive 
polarization. 

The presence of the mass could also have important applications in cos- 
mology because the fact that gravitational waves can have mass could give a 
contribution to the dark matter of the Universe. 

At the end of the paper an upper bound for these reHc gravitational waves, 
which arises from the WMAP constrains, is also released. 

2 A weak modification of general relativity 

Let us consider the action 



Equation ((T|) is a particular choice in f{R) theories of gravity [9} \W \ [TT | [T3 t 
flOl [21] in respect to the well known canonical one of General Relativity (the 
Einstein - Hilbert action [3ll[32]) which is 



Criticisms on f{R) theories of gravity arises from the fact that lots of such 
theories can be excluded by requirements of Cosmology and Solar System tests 
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|33) . But, in the case of the action ([T]), the variation from standard General 
Relativity is very weak, becauese e is a very small real parameter, thus, the 
mentioned constrains could be, in principle, satisfed ^3]. Note: General Rela- 
tivity is obtained for e = and /o = 1. 

The action ^ has been analyzed in [35 in a cosmologic context. But, 
because we will interact with gravitational waves, i.e. the linearized theory in 
vacuum. Cm ~ will be put and the pure curvature action 



(3) 



will be considered. 



3 The field equations 

Following [32l [M| (note that in this paper we work with 8ttG = 1 , c = 1 and 
h= 1), the variational principle 



6 / d^x^foR 



l+e _ 



(4) 



in a local Lorentz frame will be used. 
One gets: 



S J d''x^foR'+' ^ J d^x[S^foR'+' + foV^SR^+'] = 

= / d^x[^fo{l + e)R'R^u - i5M./oi?'+1^5^''+ (5) 
+ / d^x^{l + e)hR'cr5R,,,. 

Recalling the relation between the Christoffel coefficients and the Ricci tensor 
3 [Ml one can write 



where 



In this way, the second integral in equation ([5]) can be computed as 



(6) 
(7) 



/ d^x^{l + e)foR'g^''5R^, - / d^x^{l + e)foR'dc.X'' = 
= / d^xdc^i^il + e)hR'X"] - J d^xd^[y^{l + e)foR']X- 
Assuming that fields are equal to zero at infinity |32l [35] , one gets 



(8) 



d'^x^il + e)hR'g^''5Rf.. = - J d'^xd^^Hl + s)foR']X' 



(9) 
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Now, let us compute X". Recalling that in a local Lorentz frame it is 



V/3 9tJ.y = dpg^„ = (10) 
and using the well known definitions of the Christofell coefficients [Ml [35] it 



\g'^°'{d^,5gi3^ + d^Sg^p - dpSg^,^). 



In the same way it is 



5^1. = \9''''d,{5g,p). 



From eqs. ifTTj) and l|12p one gets 



ci'^iST^) = \d^{g^.Jgn - d^{gp^.5g-P) 



and 



(11) 

(12) 

(13) 
(14) 

(15) 
(16) 



5'^"('5r;;j = --d"{g.05g^^). 
Then, substituting in ([7]), it is 

X^^d^{g^,5gn~d^{9^,5gn. 
With this equation, equation ([9]) becomes 

/ d^x^{\ + e)foR'gt'''SR^, = 

= \ d^xdo\^){\+e)hWW{g^,Jg''-) ~ 9"(,g^,<5g'-)], 
which also gives 

/ d^x^g{\+e)hWgi'-8R^,, = 

= / d^x{g^,,d'-d^\^g(\ + e)/oi?^](5,9^''} - / d^x{g^,a'^9„[V^(l + e)hR'\bg'--}. 

(17) 

Putting this equation in the variation ([5]) one obtains 

8 J d^x^il + e)foR' ^ J d^x[^{l + e)foR'R^, - 55M./oi?'+1'5g^"+ 
+ J d^x{g^,d''d^[^{l + e)foR')] - gc.d^^d^i^il + e)foR']}Sgn 
The above variation is equal to zero for 



(18) 



(1 + e)/oi?"i?M. - 25/../oi?'+" = (Vm V. -ffM-n)(l + £)/oi?", (19) 
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which are the modified Einstein field equations. Writing down, exlplicitly, 
the Einstein tensor eqs. lfT9|) become 



(20) 

Taking the trace of the field equations (|20)) one gets 

□ (1 + e)hR' = ii^/oi?i+^ (21) 
Now, we can define the effective scalar field 

$ = (1 + e)foR' (22) 
with associated an effective potential 

dV _il-s)^ 



-foR'+'- (23) 

3q. lI2I]), a Klein - Gordc 

field is obtained 



d$ 3 

Thus, from eq. l(2T|) . a Klein - Gordon equation for the effective $ scalar 



n^ = %. (24) 
4 The linearized theory 

To study gravitational waves, the Hnearized theory has to be analyzed, with 
a little perturbation of the background, which is assumed given by a near 
Minkowskian background, i.e. a Minkowskian background plus <I> = $o (the 
Ricci scalar is assumed constant in the background) 0[19]. We also assume $o 
to be a minimum for the effective potential V: 

V ~ laS^^ m^S^, (25) 

and the constant m has mass dimension. 
Putting 

(26) 

$ = $0 + S^- 

to first order in /i^i^ and (5$, calling R^^pa , Rpu and R the linearized quan- 
tity which correspond to R^^pa , Rpv and R, the linearized field equations are 
obtained [H [H [31] : 

(27) 
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where 



(5$ 

= (28) 

<i>Q 

Then, from the second of eqs. (|27| . one can define the mass Hke 



ahm ns<P U5R^ 

Thus, as the mass is generated by variation of the Ricci scalar, we can say 
that, in a certain sense, it is generated by variation of spacetime curvature, 
re-obtaining the same result of [U [19] . The difference with the works [HI is 
that now the theory is more suitable as the modification of General Relativity is 
very weak and in agreement with requirements of Cosmology and Solar System 
tests^[33]. 

Rfii^pa and eqs. (|27p are invariants for gauge transformations [9l \T2\ fTO] 

(30) 

8^ (5$' = (5$: 



then 



can be defined, and, considering the transform for the parameter e 

□e, = ^^'h^,, (32) 
a gauge parallel to the Lorenz one of electromagnetic waves can be choosen: 

df'h^, = 0. (33) 
In this way field equations read like 

Dh^, = (34) 

Oh,n = ni^hm (35) 
Solutions of eqs. l(34| and (|35)) are plan waves [T2t[T9]: 



hf_c^ = A^i,{']f) exp{ip"xa) + c.c. (36) 
hm = a('P^) expiiq'^Xa) + c.c. (37) 



where 
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(38) 



In eqs. ((34)) and l(36| the equation and the solution for the standard waves 
of General Relativity [3ll [32] have been obtained, while eqs. l[35|) and l(37| 
are respectively the equation and the solution for the massive mode (see also 

[Smiillg]). 

The fact that the dispersion law for the modes of the massive field hm is not 
Hnear has to be emphatized. The velocity of every "ordinary" (i.e. which arises 
from General Relativity) mode hf^u is the light speed c, but the dispersion law 
(the second of eq. l(38|) ) for the modes of hm is that of a massive field which 
can be discussed like a wave-packet [9l [121 [El- Also, the group-velocity of a 
wave-packet of hm centered in "p*' is 

= (39) 

which is exactly the velocity of a massive particle with mass m and momen- 
tum 7?. 

Prom the second of eqs. (|38l) and eq. l(39)) it is simple to obtain: 



\/ up- — VP? , , , 

VG = . (40) 

Then, wanting a constant speed of the wave-packet, it has to be [51 fT2j [T9] 

= v/(1-«gV- (41) 



The relation (f4T|) is shown in fig. 1 for a value vg = 0.9. 

Now, the analisys can remain in the Lorenz gauge with trasformations of 
the type De^ = 0; this gauge gives a condition of transversality for the ordinary 
part of the field: k^^A^,, = 0, but does not give the transversality for the total 
field hfj^^. From eq. l(3T|) it is 

hf^i- ~ h^i, 'n^^f^rn- (42) 



(43) 



At this point, if being in the massless case [9l[T2l[l9], it could been put 

□e^ = 

d ff^ — —tk + h 



which gives the total transversality of the field. But in the massive case this 
is impossible. In fact, applying the Dalembertian operator to the second of eqs. 
l(43ll and using the field equations l(34|) and (|35l) it results 

□e^ = m^hm, (44) 
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Figure 1: the mass-frequency relation for a massive GW propagating with a 
speed of 0.9c : for the mass it is IHz = 10~^^eV 
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which is in contrast with the first of eqs. (|43| . In the same way it is possible 
to show that it does not exist any linear relation between the tensorial field h^i, 
and the massive field h^. Thus a gauge in wich /i^^ is purely spatial cannot be 
chosen (i.e. it cannot be put h^Q = 0, see eq. (|42|) ) . But the traceless condition 
to the field h^j,^ can be put : 

□e^ = 

(45) 

These equations imply 

= 0. (46) 
To save the conditions dfji^'^ and h^O transformations Hke 

□e^ = 

(47) 

^^,e^' = 

can be used and, taking ~jf in the z direction, a gauge in which only An, 
A22, and A12 = A21 are different to zero can be chosen. The condition h = 
gives All = —A22. Now, putting these equations in eq. I|42p . it results 

/i^^(t, z) ^A+{t- z)el+J +A''{t- z)e(,^) + /i™(i - VGz)r]^u. (48) 

The term A+(t — z)e|it^ +A^ (i — z)e^^^ describes the two standard polariza- 
tions of gravitational waves which arise from General Relativity, while the term 
hmit — VGz)r]^i, is the massive field arising from the high order theory. In other 
words, the function R"^ of the Ricci scalar generates a third massive polarization 
for gravitational waves which is not present in standard General Relativity. 



5 The primordial production of the third polar- 
ization 

Now, let us consider the primordial physical process, which gave rise to a char- 
acteristic spectrum flg^ for the relic GWs. Such physical process has been 
analyzed in different ways: respectively in refs. |251 [26] and \27\ [28] but only 
for the components of eq. (|48l) which arises from General Relativity, while in 
[29] the process has been extended to scalar-tensor gravity. Actually the pro- 
cess can be furtherly improved showing the primordial production of the third 
polarization of eq. l(48|) . 

Before starting with the analysis, it has to be emphasized that, considering 
a stochastic background of GWs, it can be characterized by a dimensionless 
spectrum [25l [26l [27l [28l [29] 



_ 1 dpg„ 

Pc din/' 



nuf) ^ (49) 
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where 

. '4 (ao) 

is the (actual) critical density energy, pc of the Universe, Hq the actual value 
of the Hubble expansion rate and dpg^ the energy density of relic GWs in the 
frequency range / to f + df. 

The existence of a relic stochastic background of GWs is a consequence 
of generals assumptions. Essentially it derives from a mixing between basic 
principles of classical theories of gravity and of quantum field theory. The strong 
variations of the gravitational field in the early universe amplifie the zero-point 
quantum oscillations and produce rehc GWs. It is well known that the detection 
of relic GWs is the only way to learn about the evolution of the very early 
universe, up to the bounds of the Planck epoch and the initial singularity [2T | [25 t 
HEl nil mi [29] . It is very important to stress the unavoidable and fundamental 
character of this mechanism. The model derives from the infiationary scenario 
for the early universe [36l [37] , which is tuned in a good way with the WM AP 
data on the Cosmic Background Radiation (CBR) (in particular exponential 
infiation and spectral index « 1 [38t [39]). Infiationary models of the early 
Universe were analysed in the early and middles 1980's (see [36] for a review 
), starting from an idea of A. Guth [37]. These are cosmological models in 
which the Universe undergoes a brief phase of a very rapid expansion in early 
times. In this context the expansion could be power-law or exponential in time. 
Infiationary models provide solutions to the horizon and fiatness problems and 
contain a mechanism which creates perturbations in all fields. Important for our 
goals is that this mechanism also provides a distinctive spectrum of relic GWs. 
The GWs perturbations arise from the uncertainty principle and the spectrum of 
rehc GWs is generated from the adiabatically-amplified zero-point fiuctuations 
[2Tl[25l[l6l[22l[28l[29]. 

Now, the calculation for a simple infiationary model will be shown for the 
third polarization of eq. (|48| , following the works of Allen [25} [26] that per- 
formed the calculation in the case of standard General Relativity and Corda, 
Capozziello and De Laurentis [221 [30] that extended the process to scalar GWs. 
Even here we have to recall that, in a recent paper [40], such a process has 
been applied to the theory arising from the action ([1]). But, in [40] a different 
point of view has been considered. In that case, using a conform analysis, the 
authors discussed such a process in respect to the two standard polarizations 
which arises from standard General Relativity. In the following the analysis is 
focused to the third massive polarization. Thus, in a certain sense, one can say 
that the present analysis is an integration of the analysis in [40j . 

It will be assumed that the universe is described by a simple cosmology in 
two stages, an infiationary De Sitter phase and a radiation dominated phase 
[m [Ml [261 [221 [281 [29]. The line element of the spacetime is given by 



ds^ = a^{T^)[~dr]^ + d-^^ + h^^{r], -^)dx^'dx'']. (51) 
In this line element, because we are considering only the third polarization. 
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the metric perturbation l(48|) reduces to 



where 

10 
r = 1 

1 ^^'^^ 
1. 

In the De Sitter phase (77 < r/i) the equation of state is P = —p = const, 
the scale factor is 0(77) = r]friQ^{2r]i — 77)"^ and the Hubble constant is given by 
H{r]) = Hds = crjo/rif. 

In the radiation dominated phase (77 > 771) the equation of state is P = p/3, 
the scale factor is a{r]) = 7//7/0 and the Hubble constant is given by H(r]) = 
cvo/v^ [2T1I251I261I2S1I3Q]. 

Expressing the scale factor in terms of comoving time defined by 

cdt — a{t)dri (54) 

one gets 

a{t) (X cxp{Hdst) (55) 
during the De Sitter phase and 

a{t) cx Vt (56) 

during the radiation dominated phase. In order to obtain a solution for the 
horizon and flatness problems it has to be [Ml l37] 

o-ivo) ^ in27 
aim) > 

The third polarization generates weak perturbations h^^iji, 'x') of the metric 
l(52l) that can be written in the form 

h^^ ^ I^,y{k)X{ri)cxp{k ■^), (57) 
in terms of the conformal time 77 where /c is a constant wavevector and 

hmiri, k , "af) ^(77) exp( k ■ "af). (58) 

By putting ^(77) = a{7])X{7]) and with the standard Hnearized calculation in 
which the connections (i.e. the Cristoffel coefflcents), the Riemann tensor, the 
Ricci tensor and the Ricci scalar curvature are found, from Friedman Hnearized 
equations it is obtained that the function Y{r/) satisfies the equation 

y" + (|fc'|2-^)y = (59) 
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where ' denotes derivative with respect to the conformal time. Cleary, this 
is the equation for a parametrically disturbed oscillator. 

The solutions of eq. l(59|) give us the solutions for the function X{t]), that can 
be expressed in terms of elementary functions simple cases of half integer Bessel 
or Hankel functions [HI [Ml [Ml [Ml [30] in both the inflationary and radiation 
dominated eras: 



For r] < rji 



for 77 > 771 



= ^[l + i/dst^"']exp-zfc(7,-7?i), (60) 



^iv) = —T~r[<^cxp —ikirj — rji) + fS expi/cfry — 771), (61) 
0(77) 

where uj = ck/a is the angular frequency of the wave (that is function of the 
time because of the Constance of fc = | A: |), a and f3 are time-indipendent con- 
stants which can be obtained demanding that both X and dX/dr/ are continuous 
at the boundary rj = iji between the inflationary and the radiation dominated 
eras of the cosmologic expansion. With this constrain it is 



" = ^+'— :^ 2:;^ 

P = ^ (63) 

In eqs. l|62l). (|63|) u = ck/a{rio) is the angular frequency that would be 
observed today. Calculations like this are referred in the literature as Bogoliubov 
coefficient methods [2T| [25l [26l [29l [30] . 

As inflation dampes out any classical or macroscopic perturbations, the min- 
imum allowed level of fluctations is that requiered by the uncertainty princi- 
ple. The solution l(60l) corresponds precisely to this De Sitter vacuum state 
[211 [Ml [Ml [Ml [30]. Then, if the period of inflation was long enough, the ob- 
servable properties of the Universe today should be the same properties of a 
Universe started in the De Sitter vacuum state. 

In the radiation dominated phase the coefficients of a are the eigenmodes 
which describe particles while the coefficients of f3 are the eigenmodes which de- 
scribe antiparticles. Thus, the number of created particles of angular frequency 
to in this phase is 



N^ = M = {^^)'- (64) 



Now, one can write an expression for the energy spectrum of the relic grav- 
itational waves background in the frequency interval (lu,lu + duj) as 

, Lo^duj , hH^ m duj hH^i H^df 



UJ 
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Eq. (|65|) can be rewritten in terms of the present day and the De Sitter 
energy-density of the Universe. The Hubble expansion rates is 

tt2 _ S-rrGpc tt2 _ S-rrGpda 

Then, defining the Planck density 

PPlanck = (66) 

the spectrum is 

Pc a In i pc a} 9 ppianck 

Some comments are needed. It has to be emphasized that the computation 
works for a very simplified model that does not include the matter dominated 
era. Including this era, the redshift has to be considered. An enlighting com- 
putation parallel to the one in [26] gives 

^^.»(/) = ^^^(l+^e,)-\ (68) 
y PPlanck 

for the waves which at the time in which the Universe was becoming matter 
dominated had a frequency higher than i?eg, the Hubble constant at that time. 
This corresponds to frequencies />(!-!- z^.qY/'^Ha, where Zeg is the redshift 
of the Universe when the matter and radiation energy density were equal. The 
redshift correction in equation (|68|) is needed because the Hubble parameter, 
which is governed by Friedman equations, should be different from the observed 
one Hq for a Universe without matter dominated era. 

At lower frequencies the spectrum is [21 } [25 t [26 t [29 l [30] 

VLg^{f)^r^. (69) 

Moreover, let us note that the results (|67|) and ([GS]) . which are not fre- 
quency dependent, cannot be applied to all the frequencies. For waves with 
frequencies less than Hq today, the energy density cannot be defined, because 
the wavelenght becomes longer than the Hubble radius. In the same way, at 
high frequencies there is a maximum frequency above which the spectrum drops 
to zero rapidly. In the above computation it has been impHcitly assumed that 
the phase transition from the infiationary to the radiation dominated epoch is 
istantaneous. In the real Universe this phase transition occurs over some finite 
time Ar, and above a frequency 

f (70) 

J max ^ /, \ A ' y 

a[to) At 

which is the redshifted rate of the transition, flg.w drops rapidly. These two 
cutoffs, at low and high frequencies, to the spectrum force the total energy 
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density of the relic gravitational waves to be finite. For GUT energy-scale 
infiation it is [2T|[25l[26l[29l[30l. 

— ^ « 10-12. (71) 

PPlanck 

6 Tuning with WMAP data 

It is well known that WMAP observations put strongly severe restrictions on 
the spectrum of relic gravitational waves. In fig. 2 the spectrum rigy^is mapped 
following [20]: the amplitude is chosen (determined by the ratio j^f^) to be as 
large as possible, consistent with the WMAP constraints on tensor perturbations. 
Nevertheless, because the spectrum falls off oc at low frequencies, this means 
that today, at LIGO- Virgo and LISA frequencies (indicate by the lines in fig. 
2) [20], it is 



ng^if )hjoo < 9 * 10-^^ (72) 

It is interesting to calculate the correspondent strain at « lOOHz, where 
interferometers like Virgo and LIGO have a maximum in sensitivity. The well 
known equation for the characteristic amplitude, adapted for the third compo- 
nent of G Ws can be used [20] : 

hracif) ^ 1.26 * I0-'\^)^hloong^{f), (73) 

obtaning [20] 

hmciWOHz) < 1.7 * 10-26. (74) 

Then, as we expect a sensitivity of the order of 10-^^ for our interferometers 
at w lOOHz, we need to gain four order of magnitude. Let us analyze smaller 
frequencies too. The sensitivity of the Virgo interferometer is of the order of 
10-21 at sa lOHz and in that case it is [20] 

hmciWHz) < 1.7 * 10^25 (75-) 

The sensitivity of the LISA interferometer will be of the order of IO-22 at 
10^"^ ~ Hz and in that case it is |20j 

hmciWOHz) < 1.7 * 10^21^ (76) 

Then, a stochastic background of reHc gravitational waves could be in prin- 
ciple detected by the LISA interferometer. 

We emphasize that the assumption that all the tensorial perturbation in the 
Universe are due to a stochastic background of GWs is quit strong, but our 
results l(72|) . l(74|) . I|75p and (|76| can be considered like upper bounds. 
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Figure 2: adapted from C. Corda - Mod. Phys. Lett. A No. 22, 16, 1167-1173 
(2007). 

The spectrum of relic sGWs in inflationary models is flat over a wide range 
of frequencies. The horizontal axis is log^Q of frequency, in Hz. The vertical 
axis is logjo V,gsw ■ The inflationary spectrum rises quickly at low frequencies 
(wave which rentered in the Hubble sphere after the Universe became matter 
dominated) and falls off above the (appropriately redshifted) frequency scale 
fniax associated with the fastest characteristic time of the phase transition at 
the end of inflation. The amplitude of the flat region depends only on the energy 
density during the inflationary stage; we have chosen the largest amplitude 
consistent with the WMAP constrains on scalar perturbations. This means 
that at LIGO and LISA frequencies, r^g.^, (/)/iioo < 9 * 10"^^ 

7 Conclusions 

It has been shown that a weak modification of general relativity produces a third 
massive polarization of gravitational waves and the primordial production of this 
polarization has been analysed adapting the adiabatically-amplified zero-point 
fluctuations process at this case. 

The presence of the mass could also have important applications in cos- 
mology because the fact that gravitational waves can have mass could give a 
contribution to the dark matter of the Universe. 
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At the end of the paper, an upper bound for these reUc gravitational waves, 
which arises from the WMAP constrains, has also been released . 
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